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Abstract. When animals are transported and pass through customs, some of them may have 
dangerous infectious diseases. Typically, due to the cost of testing, not all animals are tested: 
a reasonable selection must be made. How to test effectively whilst avoiding costly disease 
outbreaks? First, we extend a model proposed in the literature for the detection of invasive 
species to suit our purpose, and we discuss the main sources of model uncertainty, many of 
which are hard to quantify. Secondly, we explore and compare three decision methodologies 
on the problem at hand, namely, Bayesian statistics, info-gap theory and imprecise probability 
theory, all of which are designed to handle severe uncertainty. We show that, under rather 
general conditions, every info-gap solution is maximal with respect to a suitably chosen im- 
precise probability model, and that therefore, perhaps surprisingly, the set of maximal options 
can be inferred at least partly — and sometimes entirely — from an info-gap analysis. 



1. Introduction 

This paper concerns the inspection of imported herds of animals for signs of known or un- 
known major exotic infectious diseases. Imports and exports of animals represent a significant 
contribution to UK agriculture. Even though imports are subject to strict controls at the UK 
border under EU and national rules, there is a real risk of animal diseases being introduced. 
Fevre et al. [5] review the problems associated with animal movement and the spread of disease. 

We will build further on the work of Moffitt et al. [TS] , who study inspection protocols for ship- 
ping containers of invasive species, employing info-gap theory pQ to model the severely uncertain 
number of infested items. The aim of their study is to realistically take into account economical 
considerations (actual costs of testing, and of invasive species passing through customs), whilst 
also soundly handling the enormous uncertainty. 

A key feature of their, and also our, problem is that exact probabilities of the constituent 
events are very hard to come by [14] . This motivates the use of robust uncertainty models and 
decision tools, such as info-gaps [1] (i.e. robust satisficing) as in the original study, but also 
Bayesian statistics [3] and imprecise probabilities [TH], as we will do in this paper. 

Our study, using each of these decision methodologies, leads us to surmise a connection be- 
tween info-gap analysis and imprecise probability theory (r-minimax and maximality in partic- 
ular). We prove that the perceived connection is no coincidence, and we establish a rigorous 
theoretical link between the two approaches. 

The paper is organised as follows. Section [5] introduces the problem of animal inspection, 
defines the model, discusses various uncertainties involved, and derives an expression for the ex- 
pected loss under a simple binomial model for infection. Section [3] solves the inspection problem, 
first by Bayesian analysis, then using an info-gap model, and finally using an imprecise probabil- 
ity (or, robust Bayesian) model with maximality. These results are discussed in Section [4j where 
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we formally define an info-gap model based on a nested set of imprecise probability models, and 
establish the theoretical connections between info-gap, T-minimax, and maximality. Section [5] 
concludes the paper. 

2. Animal Herd Testing 

In this section, we extend a model, proposed by |I5j for the detection of invasive species, to 
suit our purpose: 

• we explicitly take specificity and sensitivity into account in order to allow for imperfect 
testing, 

• we take into account an additional cost term for terminating the herd in case an infection 
is detected, and 

• we model the occurrence of diseased animals in the herd as a binomial process, under a 
worst-case assumption of independence of infections between animals. 

2.1. Model Description. Consider a herd of n animals, of which m are tested — the problem 
is to choose m optimally. The uncertain number of diseased animals in the herd is denoted by d. 
The test has sensitivity — the probability that a diseased animal tests positive — equal to p, and 
specificity — the probability that a healthy animal tests negative — equal to q. 

Testing m animals costs c(m) utiles. If d diseased animals pass inspection undetected, we 
incur a cost of a(d) utiles. When at least one diseased animal is detected, then, typically, the 
whole herd is terminated, costing t(n) utiles. 

Following [15, p. 295, Sec. 3], in the numerical examples that follow, we take 

(1) c(m) = 1000 - 2000m + 1000m 2 

(2) a(d) = J ° lf d = ° (a = 10 000 000) 

la if d > 1 

Moffitt et al. [T5] consider n between 250 and 2 500, do not need to consider the cost of termina- 
tion (t(n) — 0), and assume perfect testing (p = q = 1). For our problem, in numerical examples 
that follow, we take 

(3) n = 250 

(4) t(n) = 400n= 100 000 

(5) p = 0.9999 

(6) q = 0.999 

so we assume that a diseased animal tests positive with probability 0.9999, and a healthy animal 
tests negative with probability 0.999. For reference, if q = 0.999, then probability that all animals 
in a healthy herd of size n — 250 test negative is q n = 0.78. 

2.2. Model Uncertainties. Obviously, many of these values are rather uncertain. The only 
values we are pretty certain of are the number of animals n in the herd, the cost of testing c(n), 
and the cost of termination t(n). 

Due to the necessity that the herd must have valid health documentation, we would expect 
that the number of infected animals d would be low. Additional inspection by veterinary officials 
is costly and depends on the inspecting official's ability to spot signs of infectious disease like 
pathological lesions and abnormal behaviour. Of course, the level of experience and competency 
will vary from official to official, but the testing procedure should be thorough enough for us to 
be confident of both a high sensitivity, p, and specificity, q. In addition to this, the government 
would prefer the most sensitive test possible (within budgetary constraints), even if specificity 
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was slightly compromised, because a rare false positive would be better for the prevention of 
disease entry than a rare false negative. Hence, we would expect p > q. Further discussion of 
this can be found in [25] . 

Of course, in general, having values for p and q as high as 0.9999 and 0.999 is unrealistic. For 
most tests, the developers have aimed at getting a high value for the specificity and sensitivity 
suffers. However, there are examples in animal disease testing where both the sensitivity and 
specificity are this high. For example, the virus antibody test for caprine arthritis-encephalitis 
claims sensitivity and specificity values of over 99.5% [T2] and near perfect sensitivity and speci- 
ficity have been estimated for the polymerase chain reaction test for parasites in fish [6]. 

Regarding the cost a of an infection passing through customs, some historical data is available. 
For example, instances of major disease outbreaks in the last couple of decades include BSE where 
public spending was over £5 billion, and the foot and mouth outbreak in 2001 which costed the 
UK government £2.6 billion [5]. These experiences show that there is great variation in the level 
of costs of exotic disease outbreaks. Due to the exceptional nature of the outbreaks, there is 
limited evidence on which to base cost assessments. Therefore, there is great uncertainty about 
what may happen in the future. 

Outbreaks of any particular exotic disease are generally rare or may never have occurred at 
all. Also, diseases change as new strains develop; consequently, the possibility of new diseases 
arriving can change rapidly. For example, until a few years ago, bluetongue was considered 
extremely unlikely in some European countries, but now outbreaks are expected every couple of 
years. 

In late 2009, an elicitation exercise was carried out with government experts to help quantify 
the average annual costs to the UK government of exotic infectious disease outbreaks and the 
uncertainty about those estimates [TT]. In that exercise, it was clear that the costs are severely 
uncertain even when the disease was known (for example, foot and mouth is an exotic infectious 
disease) . A major contributor to the uncertainty about the overall cost was the possibility of an 
outbreak of an unknown infectious disease, which could cost anywhere from £0.5 billion to £6 
billion. 

The scale and costs of an outbreak will depend on the length of time between the diseased 
animal entering circulation and the disease's presence being confirmed, and the speed and effec- 
tiveness of the government's response. The eventual costs are influenced by any public health 
implications and the effects of disease controls on other industries. The main elements of the 
costs due to control measures include: the disposal of and payments for culled animals; the 
tracing, testing and diagnosis of animals; the cleaning and disinfection of infected premises; and 
administrative costs in managing the outbreak. The size of these costs will vary according to 
the scale of the outbreak with key factors being the number of infected premises, the numbers 
of animals culled, and the duration of the outbreak. These types of factors are considered in 
greater detail in [5] and [5]. 

A serious study of how all uncertainties involved could be taken into account in the model 
would of course be extremely interesting, but is beyond the goal of this paper. Instead, in this 
initial study, following |15) and many others, for now we will focus on the main uncertainty, that 
is, the number of diseased animals d, and simply assume reasonable values for the remaining 
parameters. Such simplistic model helps to illustrate the methodological differences and to 
motivate the theory. 

2.3. Expected Loss. First, we derive the expected loss, in case all parameters of the problem 
are perfectly known, including the number of diseased animals d. Clearly, conditional on d, the 
expected loss is: 



(7) 



L(m,d,p,q,c,a,t) = c(m) + t(n) Pr(T|d) + a(d) Pr(T c |d) 
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Figure 1. Loss as a function of the number of diseased animals for to = 10 and 
to = 20. 



where T denotes termination of the herd, that is, the event that at least one diseased animal is 
detected, and T c denotes its complement, that is, the event that the herd passes inspection. 

Let us deduce Pr(T c |d). First, if the test group of size to is sampled randomly and with- 
out replacement, then the probability of exactly z diseased animals in the test group follows a 
hypergeometric distribution: 

/d\ tn—d\ 

(8) Pr(z|d) = ^f-. 

Next, we calculate the probability of non-termination given z diseased animals in the test 
group, that is Pr(T c |d, z). If d = 0, then the probability of non-termination is the probability 
of all healthy animals in the sample testing negative, so Pr(T c |0, z) — q m . If d > 1, then given 
z diseased animals in the sample, non-termination occurs when none of the z diseased animals 
tests positive and all of the m — z healthy animals test negative. Hence, in all cases, 



(9) 



Pr(T c |d,z) = (l-p)V 



By the law of total probability, 

Pr(T c |d) 

(10) 



d 

£ 

z=0 



Pr(T c |d,z)Pr(z|d) 



(f)( 



n—d\ 



)ji ,,)-<r > - ,> '" 



z=0 



Now we have all the ingredients to calculate the total expected loss if we choose to test to out 
of n animals: 

(11) L(m,d,p,q,c,a,t) = c(m)+t(n) + (a(d) - t(n)) Pr(T c |d) 

or, if a'(n,d) = a(d) — t(n) denotes the termination adjusted cost of apocalypse, 

(12) = c(to) +t(n) +a'(n,d)Pi(T c \d) 

where Pr(T c |c?) is given by Eq. (fit))) . Figure []] depicts the expected loss for a few typical cases. 
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2.4. A Binomial Model for Infection. Moffitt et al. [TS] consider an info-gap model directly 
over the number of diseased animals d, which leads to a rather tricky optimisation problem. 
Instead, we will consider the (highly uncertain) probability r that an animal is infected, and 
derive the expected loss as a function of r. Although we do not explore this topic further in this 
paper, this also paves the way to modelling spatial dependencies between infections in the herd, 
leading to more optimal testing strategies. 

So, assume that each animal has a probability r of being infected; for simplicity, for now, we 
assume that one animal being diseased does not affect another animal being diseased. Obviously, 
this will generally not be satisfied, and more realistically, we would expect a positive correlation, 
resulting in diseased animals being clustered together in the herd. Assuming independence 
essentially amounts to a worst case study: at the other extreme end, if one diseased animal 
would immediately infect the whole herd, then it would be sufficient to test only a single animal, 
as d — and d — n would be the only two possibilities. 

Under the worst case assumption of independence, the probability of having d out of n animals 
infected is: 

(13) Pr(d|r) = (^)r d (l - r) n - d 
The expected loss is: 

a 

(14) E(L(m, -,p, q, c, a, t)\r) = ^ L(m, d,p, q, c, a, t) Pr(d|r) 

From now onwards, we will simply write L(m\r) instead of E{L(m, -,p, q, c, a, t)\r) in order to 
simplify notation. Figure [2] depicts L(m\r) as a function of m for a few typical situations. 

3. Decision Analysis 

In this section, we explore and compare three decision methodologies, designed for severe 
uncertainty, on the problem at hand. In particular, 

• we perform a Bayesian analysis, 

• we accommodate the info-gap approach suggested by [TS] to our extended model, 

• we investigate possible ways of constructing sets of probabilities (i.e. imprecise proba- 
bility models) which are in some sense equivalent to the proposed info-gap model, and 

• we compare the decisions that these various models lead to. 

3.1. Bayesian Analysis: Minimising Expected Loss. In this approach, we model the un- 
certainty about d probabilistically and choose the m that minimises expected loss. Of course, 
this is not the only strategy open to us when we use a Bayesian formulation, but it is the most 
common. 

Using the same model for d as in Eq. (|13p . we can write the probability mass function for d 
to be 

(15) Pr(d) = / Pr (d\r)p(r)dr, 



where p(r) is a probability density function chosen to characterise our beliefs about r. For 
computational convenience and because the density is sufficiently flexible, we choose to model 
r ~ Beta(a, /3), where a and /3 are chosen so that the Beta density matches our beliefs. Using 
this density for r, 

i\ Be(a + d, (3 + n — d) 



(16) Pr(d|Q,/3)= 



Be(a,/3) 
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Figure 2. Expected loss L(m\r) as a function of the test group size m, for 
r = 0.00010, r = 0.00025, r = 0.00050, and r = 0.00100, from bottom to top. 
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Table 1 . Choices for parameters of the Beta distribution over r. For informa- 
tion, we also state the standard deviation a, and the a and /J parameters of the 
canonical parametrization. 



where Be(-, ■) denotes the usual Beta function. In this analysis, we consider various choices of 
parameters for the Beta distribution, as listed in Table [TJ Following Walley [22 , t denotes the 
expectation of the prior, and s denotes a + (3, so a — st and /3 = s(l — t). In our analyses, 
we have varied the expectation of r by varying t, and we have chosen s such that the standard 
deviation is 0.001 throughout. We have chosen to investigate this set of distributions because 
they cover a range of reasonable beliefs that we may have about r when we are dealing with rare 
diseases. 

For each Beta distribution, we can calculate the expected loss for each choice of m using 



(17) 



n 

E{L\a,0) = Pr(d\a, f3)L(m, d,p, q, c, a, t). 

d=0 
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Figure 3. Expected loss plotted against m for each choice of parameters. Lower 
curves correspond to lower values of t. 
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Figure 4. Probability of exceeding loss threshold L c for m = 10 and r ~ 
Beta(a, (3) with parameters corresponding to the worst case considered (bottom 
row of Table [T]), as a function of the critical cost L c . 



Figure [3] is a plot of the expected loss against number of tests for the three chosen distributions. 
The minimum expected loss occurs for m = m* , given in Table [T] 

Also, by taking a probabilistic viewpoint, we can derive a distribution for the possible losses 
given a choice of m and p(r). If we choose m = 10 and the worst case distribution for r (bottom 
row of Tabled] where there is a greater than 30% chance of there being diseased animals), we can 
find probabilities for the losses exceeding different values. For the thresholds in Figure^] as soon 
as we exceed the cost of termination and testing, we are essentially calculating the probability 
of disease outbreak. The jumps in this plot correspond precisely to the cost of testing (81 000), 
and the cost of testing plus termination (181 000). 

It is worth noting that these probabilities are particularly sensitive to the choices we make for 
the parameters of the beta distribution. However, it is straightforward to assess the influence 
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Table 2. Probability of exceeding loss threshold L c = 181 000 for m = 10 and 
r ~ Beta(a, /3) for a wide range of parameter choices. Values are shown for 
L c = 182 000, as the probability is discontinuous at L c = 181 000 (see Figure E]). 



the parameters by tabulating the probability of L c exceeding the cost of testing and termination 
(181 000 when m = 10) as shown in Table EJ 

3.2. Info-Gap Analysis: Maximising Robustness. Another approach to solve our decision 
problem, under severe uncertainty about the exact probability r of a single animal being infected, 
is to select that decision which meets a given performance criterion, L c , under the largest possible 
range of r. Given that we have almost no information about r, except that it assumes a very small 
but otherwise unknown value, this simple model seems to suffice for our purpose. Obviously, 
one could define many other more refined info-gap models — and our choice of model is just one 
example among many. For a much more detailed account, see [jQ. 

Specifically, for a given value of L c , the largest possible range [0, h] of r for which we meet 
our performance criterion is characterised by 



(18) 



him, L c ) = max < 

h>0 



h: 



max L(m\r) 

r£[0,h] 
r<l 

M(m,h) 



< U 



The value h(m, L c ), as a function of L c , is called the robustness curve: it tells us how uncertain 
about r we can be for our decision m still to meet a given level of performance L c . 

A quick Poisson approximation reveals that as long as exp(— nh) is sufficiently close to 1 (and 
this holds for sufficiently small values of nh) the inner maximum over r 6 [0, h] is achieved at 
r = h (also see Figure [2] the cost increases as r increases), so 

(19) M(m,h) = L{m\h) 

Obviously, M(m, h) increases as the horizon of uncertainty h increases, whence h(m, L c ) as a 
function of L c is simply the inverse of M (m, h) as a function of h. In other words, plotting 
M{m, h) as a function of h for different values of m effectively gives us the robustness curves. 
Figure [5] depicts them. 

The choices of m which maximise robustness, for various values of the critical cost L c , are 
tabulated in Table[3] For example, at an expected cost L(m\r) of at most L c — 3 000 000, we can 
safeguard against any probability of infection r £ [0, 0.001 479], by testing 10 animals in the herd. 
For comparison, the Bayesian expected cost E(L\a,j3) for t = 0.0016 is almost exactly equal to 
3 000 000 (last row of Table[T|) when m — 10, which is in agreement with the info-gap analysis. Of 
course, it is to be noted that the actual inputs into each decision model arc different; nevertheless, 
both are assuming extreme uncertainty of a similar order if we interpret the Bayesian analysis 
as a worst case analysis — which we can easily do here due to the simplicity of the model. Thus 
the agreement is not all that surprising. 
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Figure 5. Robustness curves h(m,L c ) as a function L c for test group sizes 
m = 1 (solid), m = 15 (dashed), and m — 30 (dotted). 



ic/10 6 


m* 


W 3 h(m*,L c ) 


0.5 


2 


0.207 


1.0 


4 


0.426 


1.5 


5 


0.661 


2.0 


6 


0.912 


2.5 


8 


1.184 


3.0 


10 


1.479 


3.5 


11 


1.803 


4.0 


13 


2.163 



Table 3. Info-gap choice of m, and corresponding horizon of uncertainty, for 
various values of the critical cost L c . 



3.3. Imprecise Probability Analysis: Maximality over a Partial Ordering. There are 
several ways one might go about constructing an imprecise probability model for our problem. 
As we have just seen, the info-gap approach hinges on the idea of satisficing: we may start out 
with a level of minimum performance that we hope to achieve, and the analysis tells us how 
much uncertainty we can tolerate, at this price. One might also interpret it conversely: for a 
given level of uncertainty, the analysis tells us how much we might potentially pay, if it comes 
to the worst. 

Typical decision models for imprecise probabilities studied in the literature do not relate to 
satisficing, yet, they do incorporate an idea similar to the info-gap horizon of uncertainty: the 
imprecision of our model. Concretely, consider the set M.h of all probability densities over r that 
are zero outside [0,/i]Q We say that a choice m dominates a choice ml, and we write m y m! 
whenever the expected loss under m is strictly less than the expected loss under m' over all 
densities p in A4h, that is, whenever 

p oo poo 

(20) / L(m\r)p(r)dr + e < / L(m'\r)p(r)dr 

Jo Jo 



The adventurous reader may take all finitely additive probability measures [i on [0, +00] with fi([0, h]) = 1. 
We do without this complication: because all functions involved are continuous, those additional measures make 
no difference. 
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Table 4. Result of Eq. (|2"2"|) (divided by a factor 10 3 for everything to fit in the 
table). A positive value means that the corresponding choice of m is optimal 
for the given horizon of uncertainty h. 



for all probability densities p in Mh and some e > 0. 

Because L(m\r) is continuous in r for every m, this happens if and only if 

(21) min \L(m'\r) - L(m\r)} > 

re[Q,h] 

Note that the min rg r w operator can be thought of as a lower expectation operator, or lower 
prevision P h — we will come back to this in Section |U 

One can easily prove that >- is a partial order, whence, a sensible way to choose m is to pick 
one which is not dominated by any other option, or in other words, which is maximal. The idea 
of choosing undominated options goes back at least to Condorcet [H pp. lvj-lxix, 4. e Exemple]; 
also see [T7J p. 55, Eq. (1)], [2U Sections 3.7-3.9], and [T^] for further discussion. Maximality 
has also been used in robust Bayesian models [131 §10.4], under slightly different terminology. 

Given our partial order, one can easily show that an option m is maximal if and only if 

(22) min max \L(m'\r) — L(m\r)] > 

ro'e{0,i,...,n}i-e[o 1 /j] 

The inner maximum is almost always achieved at either r = or r — h, simplifying practical 
calculations substantially. Table E] depicts these values for all choices of m, and varying values 
of h. For ease of comparison with the info-gap solution, we have chosen the same values of h as 
those listed in Table [3] 

4. Discussion 

It is already well known that robust Bayesian models and imprecise probability models are, 
for the most part, mathematically equivalent [2lJ §5.9, pp. 253-258]. Therefore, in the following, 
we will focus on info-gap and maximality — the latter being used for both imprecise and robust 
Bayesian models. 
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Interestingly, in our example, info-gap and maximality give essentially the same result, with 
maximality refining the picture slightly: for a given horizon of uncertainty h, the maximal 
solutions are {1, . . . ,m*}, where m* is the info-gap solution. The most notable result is that 
all info-gap solutions are maximal. Is this a coincidence? Formulating info-gap theory in terms 
of lower previsions, we show that this holds for arbitrary info-gap models and arbitrary lower 
previsions, subject to the mild and usually satisfied conditions of Theorem [T] (for T-maximin) 
and Theorem [2] (for maximality). 

4.1. Info-Gaps for Imprecise Probabilities. Let u 6 be an uncertain parameter of interest — 
can be an arbitrary set. We must select a decision d from a finite set D. The loss function 

L(d,uj) represents the loss (in utiles) if we choose d and co obtains. 

Info-gap theory starts out with a family of nested sets Uh of 51, where h is a non-negative 

parameter called the horizon of uncertainty and [/j,C Uh> whenever h < h' . In our example, Uh 

was simply [0, h]. Following that example, we saw that a very natural way to model these nested 

sets Uh in terms of sets of probabilities goes by way of a vacuous model Mh, that is, the set of 

all probability densities that are zero outside Uh- 

If we denote the upper expectation induced by Mh by Ph (i.e. the pointwise lowest upper 

bound for the set of expectation operators associated to Mh), then, formally, we define the 

info-gap solution D*(L C ) C D at satisficing level L c as: 



Note that D*(L C ) will usually be a singleton (or, the empty set). 

Also note that the first equation may not have a solution: this happens when Po(L(d, •)) > L c , 
that is, when d is infeasible even if we are as certain as can be (h = 0). 

Now, from the point of view of imprecise probability, there is no compelling reason to restrict 
ourselves to vacuous models. In fact, we can allow Mh to be any set of probability densities 
on Q — we already picked a more general set in our Bayesian analysis in Section 13.11 — under one 
restriction: a close inspection of the theory reveals that a crucial property that the info-gap model 
relies on is that the worst case cost, Ph(L(d, •)) is non-decreasing as the horizon of uncertainty 
h increases. Whence, we logically impose that Mh Mh' whenever h < h' . 

So, instead of starting out from a family of nested subsets Uh of f2, we start out from a family 
of nested sets Mh of probability densities on Q. For instance, in our example, the uncertainty 
was over the values of r, so Mh would be some arbitrary set of probability distributions for r. 
In the Bayesian analysis, we restricted Mh to the Beta family, for computational convenience. 

One can of course interpret this again as an info-gap model, where the uncertain parameter 
is now the probability density over fi — also see [H pp. 1062-1063] for an informal discussion of 
this approach. The imprecise Dirichlet model [22^ is an example of such family (with h = 1/s). 
For another example, see [7j for a discussion of nested sets of p-boxes and the resulting info-gap 
analysis. 

4.2. Main Result. The next result links the info-gap solution to the so-called r-minima>0 
solution. See [2 p. 1061, Fig. 14] for an informal discussion of a very similar equivalence between 
info-gap and minimax. Interestingly, [18l p. 4, Table 1] constructs a maximin model that is 
fully equivalent to an info-gap model. The result below is quite different, as we change neither 
variables nor loss function. 



T-minimax minimises the upper expectation of the loss 1201 110| . 



(23) 
(24) 



h{d,L c ) :=max{/i: P h (L{d,-)) < L c ) 
D*(L C ) := argmax h(d, L c ) 
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Effectively, we show that, if certain fairly mild conditions are satisfied, the info-gap solution 
coincides with the r-minimax solution — parametrized over the horizon of uncertainty h — of ex- 
actly the same problem. In fact, such approach has already been used as a technique to solve 
info-gap problems (see for instance [16l PP- 1690-1691]). Below, we identify sufficient conditions 
for this to work, and provide two counterexamples in cases where one of these conditions fails. 

For convenience, we denote the r-minimax loss at horizon h € R + by L*(/i)2 

(25) L*(h) :=min P h (L(d,-)). 

deD 

Note that L* is non-decreasing as a function of h, because each Ph(L(d, ■)) is. We will be 
interested in the right derivative of L* at a point h: 

h'^h h! — a 

h'>h 

Because L* is non-decreasing, d+L*{h) is non-negative and well defined (possibly +oo) for every 
heR+. 

Theorem 1. Let h G R + and L c € R. The info-gap solution D*(L C ) coincides with Y -minimax 
solution with respect to Ph, that is, 

(27) D*(L C ) =&rgmm P h (L(d,-)), 

deD 

whenever the following conditions are satisfied: 

(28) d+L*{h) > 0, and 

(29) L*(h)=L c 

Proof. By definition, d* € D*(L C ) whenever, for all deD, 

(30) h(d*,L c )>h(d,L c ) 
By definition of h(d, L c ), this is equivalent to saying that 

(31) {h':P h >(L(d*,-))<L c }DU deD {h':P h >(L(d,-))<L c } 
Rewriting the above expression, we have, equivalently, 

(32) {h': P h ,(L(d*,-)) < L c ) D \ ti: mmP h ,{L(d, ■)) < L c \ 
k ' \ deD J 

or, by Eq. (g5]), 

(33) {h':P h ,(L(d*,-)) < L c ) D {h': L*{h') < L c } 

By Eq. and the non-decreasingness of L* , it is easily seen that L*{h') > L*(h) for all h! > h. 
Moreover, by Eq. (|2"9"|) . L c = L*(h). Concluding, the set on the right hand side is a fancy way of 
writing [0, h]. Therefore, the above is equivalent to 

(34) P h {L{d*,-))<L c 

Once more by Eqs. (|2"9")l and (|2~5|) . this is equivalent to saying that d* is a r-minimax solution 
with respect to Ph- □ 



3 We set R+ := {x G M: x > 0}. 
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Interestingly, for given L c such that 

(35) mm P {L(d, •)) < L c < minP 00 (L(d, •)) 

d£D d£D 

it holds that Eq. (|29[) has a unique solution for h > whenever L* is strictly increasing and 
continuous in h. This means that we are effectively free to choose L c under the additional 
assumption of continuity. 

To see why we are not free to choose L c when continuity is not satisfied — and violate Eq. (j29|) — 
imagine for instance that: 

(36) *<*<*.•» * h K *\ *<*<*.•»■ 

Then, for L c — 3, we have that D*(3) = {di,^} because h(d,3) = 1 for both d\ and d%, yet 
obviously d\ is T-minimax (it could even be uniformly dominated by d 2 ). Effectively, this is 
simply a technical limitation of the info-gap model, as any reasonable person would probably 
agree with the T-minimax solution. 

What happens if Eq. ((28)) is violated? Imagine for instance that: 






(37) P h (L(d l ,-)) = i 1 . f ~ P h (L(d2,-)) = 



With this choice, L* is continuous, but not strictly increasing. For h = 1, we have that 
argmin de£) P h (L(d, •)) = {di,d 2 } because P/,(£(di, •)) = P h (L(d 2 ,-)) = for /i = 1. Yet, 
for L c = = we have that h(d\, L c ) = 1 < hid^, L c ) = 2, so only d 2 is optimal according 

to the info-gap criterion. This example uncovers a technical limitation of the r-maximin model, 
as for this case, any reasonable person would probably agree with the more robust info-gap 
solution. 

Now, it is well known that every T-minimax solution is also maximal (see for instance |19]). 
whence, we conclude: 

Theorem 2. Suppose d+L*(h') > for all h! E [0,h]. Then, for all h' € [Q,h], every info-gap 
decision d* € D*(L*(h')) is maximal with respect to Ph- 

(38) [J D*(L*(h'))C{deD:(WeD)(P h (L(d',-)-L(d,-))>0)} 

h'e[0,h] 

Proof. Use the preceding theorem, and note that every T-minimax with respect to P_ h , is maximal 
with respect to P_ h , provided that h' G [0, h]. □ 

Again, if in addition L* is continuous on [0, h], then the range for L*{h') in the above theorem 
is simply an interval: 



(39) {L*(h'):h'e[0,h]} = 



mmP (L(d,-)),minP h (L{d,-)) 

d£D d£D 



Summarising, Theorem [T] provides sufficient conditions for the info-gap solution, for fixed 
values of L c and h, to be equal to the T-minimax solution: proponents of cither approach are 
'observationally equivalent' [3J Sec. 7]. 

Theorem [5] shows that a full fledged info-gap analysis, varying the horizon of uncertainty 
along an interval, yields an elegant approach to capture maximal solutions. In our example, we 
actually find all maximal options — in general this may not be the case. Still, it shows that an 
info-gap analysis can be of value even if maximality is the final goal: 

• an info-gap analysis might give a rough idea of the size of the maximal set (in particular, 
it provides a lower bound for it), 
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• the analysis can be an appealing way to represent the maximal solution graphically (as 
in Figure [5]), and 

• as robustness curves show the trade-off between uncertainty and cost, they are also 
obviously useful in the process of elicitation. 

5. Conclusion 

We constructed a simple model for inspecting animal herds for dangerous exotic infections, 
building further on the work of Moffitt et al. |15j . We solved the problem using three popular 
decision methodologies suited for dealing with severe uncertainty: Bayesian analysis, info-gap 
analysis, and imprecise probability theory (maximality and r-minimax). We found that, in this 
example, the solutions of the info-gap and imprecise models essentially coincide, although the 
way they arrive at it is very different. 

We explored the theoretical link between info-gap theory, r-minimax, and maximality. We 
established that, under rather general conditions, every info-gap solution is maximal. Therefore, 
the set of maximal options can be inferred at least partly, and sometimes wholly, from an info- 
gap analysis. Consequently, robustness curves also make sense in an imprecise probability (or, 
robust Bayesian) context, for exploring maximal options, and for elicitation, when studying the 
trade-off between uncertainty and cost that is often of interest to decision makers. 
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